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ABSTRACT 

For purposes of stellar population analysis, emission corrections for Balmer series 
indices on the Lick index system in Sloan Digital Sky Survey (SDSS) stacked quiescent 
galaxy spectra are derived, along with corrections for continuum shape and gross stellar 
content, as a function of the Mg b Lick index strength. These corrections are obtained 
by comparing the observed Lick index measurements of the SDSS with new observed 
measurements of 13 Virgo Cluster galaxies, and checked with model grids. From the Ha 
Mg b diagram a linear correction for the observed measurement is constructed using best 
fit trend lines. Corrections for H/3, H7 and H<5 are constructed using stellar population 
models to predict continuum shape changes as a function of Mg b and Balmer series 
emission intensities typical of HII regions. The corrections themselves are fairly secure, 
but the interpretation for H<5 and H7 indices is complicated by the fact that the H<5 and 
H7 indices are sensitive to elemental abundances other than hydrogen. 

Subject headings: galaxies: abundances — galaxies: clusters: individual: (Virgo cluster) 
— methods: analytical — surveys 



1. Introduction 



One of the more important factors in determining the age of a stellar 



ellipti cal galaxy is being able to measure the hydrogen Balmer series IWorthev 



p opula tion such as an 
|l994h and lO'Conneh 



(|1976l ). The first two Balmer lines Ha and H/3 are primarily used to deter mine the age of a s tellar 
population becau se of their relative insensitivity to changes in metallicity (j Serven et al.l (|2005l ) and 
Korn et al.1 (|2005l )) making a measurement of these two lines a more reliable estimato r of age then 



the h igher o rder Balmer lines H7 and H<5, which are sensitive to changes in metallicity IServen et al 



(|2005l ) and lKorn et al.1 (|2005l ) making a knowledge of the relative abundances of elements in a given 
spectra necessary in order to get a consistent agreement of ages from all the Balmer lines. 

Another important factor in determining the age of a galaxy is hydrogen emission, which 
can give the impression of a much older galaxy by filling in and weakening the diagnostic Balmer 
absorption features. Traditionally, an emission correction is used that is based on the Oil or 
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OIII emission lines (jSchiavon et ajj[2006j ; lGonzalea ll993) in which the difference in the measured 
equivalent width(EW) of a hydrogen line such as Ha or H/3 due to emission is some fraction of the 
EW of Oil or OIII. However, the correlation between hydrogen and oxygen emission strengths is 
weak, and indeed there is no astrophysical reason why the two should be tightly correlated. On 
the other hand, being able to use hydrogen recombination lin es, which have ratios that are nearly 
fixed with respect to themselves (lOsterbrock fc Ferlandl 120061 ) should give a more reliable result. 



Toward that end, we compare with 13 high-quality spectra of Virgo cluster early-type galaxies 
(Serven Jk W orthey, in preparation) with a sample of SDSS spectra provided by G. J. Graves 
(|Gravesl 1200^ . 



The Virgo Cluster elliptical galaxies consist of 13 red galaxies with colors 0.75 <B— V< 0.97. 
They also have a high S/N from S/N=150 to S/N=450, are well fluxed using standard star flux 
calibrations, and have tight control over instrument resolution and velocity dispersion. These 
spectra are long slit spectra obtained using the T2KB chip on the Cassegrain spectrograph on 
the Kitt Peak Mayall 4 meter telescope (2006 Jan 31 - Feb 5) and were chosen to cover velocity 
dispersions from 50 to 350 km s -1 with a wavelength range from 3200 A to 7500 A in two wavelength 
swaths. 



The iGravesI (120071 ) sample was constructed by Graves from spectra of 22,501 SDSS galaxies 



that fall into the redshift range 0.06< z <0.08 and who ' s colo rs meet the criteria (g-r) 



o.i 



> 



0.25( M r °- 1 - 5 log h ) + 0.42, with h = 0.70 lYan et al.1 (j2006l ) placing them firmly within the 
red sequence. These spectra are then further divided into those with Ha and Oil emission lines 
and those without. Those with emission lines were said to be "LINER-like" and those without 
"quiescent", where LINER stands for Low Ionization Nuclear Emission-line Region. 

These quiescent galaxies, which are the object of this paper, were then reduced to a sample 
of 2000 by constructing a Gaussian distribution centered about EW Oil = A with <J\oii] = 1-56 
A to match the width of the distribution. This was then truncated at (T\oii] to exclude outliers 
that may contain emission. From this distribution 2000 quiescent galaxies where randomly selected 
and split into 6 velocity bins such that they contain roughly equal numbers of galaxies. The bins 
ranges are 70< a <120 km s _1 , 120< a <145 km s -1 , 145< a <165 km s -1 , 165< a <190 km 
s _1 , 190< a <220 km s _1 , and 220< a <300 km s _1 . In order to compare the individual bins, 
the individual spectra within each bin were smoothed to a a = 300 km s -1 and coadded to make 
six composite spectra. Along with these spectra the S/N at each resolution element was computed 
to produce an error spectrum for each of the six composite spectra. Factors that may contribute 
to the error estimate are age, individual abundances, and emission signal. The error spectrum is 
likely to be dominated almost entirely by measurement uncertainty given the poor resolution of the 
individual SDSS spectra. 

Measurement from these composite spectra and their errors are what is represented as the six 
SDSS data points in this paper's fi gures. For mo re on the binning of the SDSS spectra and more 
on the SDSS spectra in general see I Graves! (120071 ). 
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The central aim of this paper is to measure the strength of any residual hydrogen emission is 
these quiescent spectra. Noting the apparent residual Ha emission of the SDSS spectra as seen in 
the first panel of Figure 1 we construct simple formulae for the correction of the relative intensities 
of the rest of the Balmer series for the Sloan Digital Sky Survey (SDSS) as a function of Mg b 
that includes continuum slope effects and the intrinsic decrement values for the Balmer lines. The 
corrections should find future use with integrated-light models to better predict stellar populations 
parameters, especially the mean metallicities and ages of galaxies. 

Our analysis methods and results are set out in the following section with a discussion and 
summary in § 3. 



2. Analysis and Results 

The first panel of Figure 1 shows that there is a discrepancy in the Ha measurements be- 
tween the SDSS an d Virgo spectra for this work use measurements of the Ha index as defined by 



Cohen et al.1 (|1998l ). Plausibly, this is either due to hydrogen emission in the SDSS spectra or to 
nonspectrophotometric wavelength-dependent fluxing errors that would propagate into Lick index 
measurements depending on how the various passbands lay across the various spectral "lumps." 
We infer that the bulk of the discrepancy is most likely due to hydrogen emission in the SDSS 
spectra because the continuum-shape differences between the data sets is too small to account for 
much of the discrepancy, as we now show. 

To determine if differences in the continuum of the two data sets could explain this discrepancy, 
the ratio of similar spectra (similar as regards velocity dispersion) from the two data sets was 
computed and a continuum fit was found for this ratio using the "continuum" routine in IRAF 0. 
Lick style indices were measured from the fit as one would measure spectra. The delta-index thus 
discovered was found to be small: about 0.06A of the measured values for Ha, insufficient to explain 
the SDSS- Virgo discrepancy. It should be noted that the effects of the continuum differences could 
be more severe if differences in the continuum shapes exist that are of the same wavelength span 
as the features of interest. Most would likely lie with the SDSS spectra, since our fluxing of the 
Virgo data was careful, but, unfortunately we don't have the data to check the significance of these 
differences. 

We character ized the SDSS- Virgo Ha - Mg b trend by best fit lines calculated using fitexy.f 



(jPress et al.lll992l ) . This is a fortran program for finding the best fit line for data with errors in both 
the x and y coordinate. This minimizes the distance of each point form the line while taking into 
account weighting by the precision of the individual measurements in both the x and y coordinate. 
The choice to characterize the trends in Ha as a function of Mg b was for two reasons. One was 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under Cooperative agreement with the National Science Foundation. 
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the tight correlation between Mg b and a and the other was for when using models as described 
later Mg b was just practical. Those line fits and the root mean squares (RMS) of the distances 
of the points from their fit lines. Also shown in Table 1 is the form of the correction term (j^) 
determined from the following derivation. 



Table 1 



Data set 


Line Fit in A 




RMS of Fit in A 


Virgo 


3.0132 - 0.3768*Mg 


b 


0.066 


SDSS 


1.3476 - 0.0532*Mg 


b 


0.020 


Correction Term 


1.6656 - 0.3236*Mg b = 


_ ja 
Fc,a 


0.069 



Table 1: Shown in this table are the best fit Ha vs. Mg b lines for the Virgo and Sloan data sets 
and their associated RMS values (lines 1 and 2). Line 3 is the difference between these two best fit 
lines, which represents the Ha correction term along with it's RMS value. 



Using the linear correction term for the hydrogen emission in Ha a correction for the subsequent 
Balmer lines was constructed under the assumption that the entire shift is due to Balmer emission 
fill-in. To determine the form of the correction term we start with the definition of equivalent width 
(EW; see Eq. 1). In Eq. 1, Ai and A2 are defined as the blue and red wavelength bounds of the 
ind ex passband an d Fj an d F c are the average index flux and the pseudo continuum flux respectively 



sec 



Worthev et all ()1994h . 



EW = (1 - -±)d\ = AA(1 - (1) 

J Al c c 

With more generality, and including a term for the flux due to an emission feature, 

EW = [ X \l - Fx ^ Fj )d\ = AA(1 - ^ - §) (2) 



Al 



F c ' v F c F c 



where Fj is the flux of the emission feature and Fx is the flux of the stellar light and AA = 
A2 — Ai. One should note, however, that it is the emission line's power that is constant. If we call 
the line's power j then the average emission line flux is defined by 

r\2 

j = FjdX = AXFj (3) 

Thus, the average stellar flux inside the continuum band is F{ = F\ and the correction term 
for the equivalent width is AXFj/F c = j/F c . To extend to Balmer lines other than Ha, we exploit 
the fact that the decrements jp/ja, jj/ja, and js/ja are known, and relatively constant. 
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For example, if j a is known, then extending from an Ha index to an H/3 index is accomplished 
by adding the correction term 



<c,a J/3 



J3 3c 



(4) 



Having determined the Ha correction from best line fits (see Table 1) all that is left to do is 
determine the conversion factors and 



For this paper a version of the lWorthevI (119941') andlTrager et al.l (|1998l )) models were used that 
use a grid of synthetic spectra in the optical iLee et al.l (J2009J) in order to investigate the effects of 
changing the detailed elemental composition on an integrated spectrum was used to create synthetic 
spectra at a variety of ages and metalliciti e s for single-burst stellar populations. The underlying 
isochrones for this paper were the IWorthevI (| 19941 ) isochrones, because they allow us "manual" HB 
morphology control. However, there are certain caveats to using these isochrones. specifically the 
models are a bit crude by today's standards and the ages are about 2 Gyr too old, so that 17 Gyr 
should really be interpreted as 15 Gyr. Other isochrone sets were used to check the results. 

For this paper new index fitting functions were generated. The data sources include a variant 
of the original lick collection of stellar spectra (jWorthev et al.lll994l ) in which the wavelength scale 
of each observation has been refined via cross-correlation, as well as the MILES spectral library 
(ISanc hcz-Bl azquez et al.ll2006r ) with some zero point corrections, and the Coude Feed Library (CFL) 
of lValdes et al.1 (120041 ). The CFL was used as the fiducial set, in the sense that any zero point shifts 
between libraries were corrected to agree with the CFL case. The MILES and CFL spectra were 
smoothed to a common Gaussian smoothing corresponding to 200 km s . The rectified-Lick 
spectra were measured and then a linear transformation was applied to put it on the fiducial 
system. 

Multivariate polynomial fitting was done in five overlapping temperature swaths as a function 
°f Qp// = 5 40/Tp .f f , log g, and [Fe/H]. The fits were combined into a lookup table for final use. As 
in IWortheyl (| 19941 ). an index was looked up for each "star" in the isochrone and decomposed into 



"index" and " continuum" fluxes, which added, then re-formed into an index representing the final, 
integrated value after the summation. This gives us empirical synthetic spectra when variations 
in chemical composition are needed. The grid of synthetic spectra is complete enough to predict 
nearly arbitrary composition. 

These models were then used in order to determine the conversion factors for the continuum 
levels. By first taking the ratio of the continuum levels near Ha and the rest of the Balmer lines 
as measured from these models and plotting them against Mg b. Then, fitting a least squares line 
to the data. This line then gives the change of the relative continuum levels as a function of Mg b. 



The last piece of the puzzle is that the native relative intensities of the Balmer lines need to 
be accounted for. This is taken into account by the relative Balmer line intensities as calculated 



by lOsterbrock Ferlandl (|2006l ). We adopt case B, 10000 K conditions, because they are near the 
middle of the range for star formation regions (j a / ja = 2.85), an d they are not too drastically 
different than LINER type spec tra in/is = 3.27 from Graves! ([20071 ). The final correction formulae 



were applied to H/3 defined in IWorthev et all (jl994l )7 HjA, H/yjr, H5a, and H5f as defined in 
(jWorthev fc Ottavianilll997l ) and are found in Table 2. 



Tabic 2 



B aimer 


Correction 


Ha Correction 




Continuum Correction 


Decrement 


Correction term 


Index 


Term in A 


in A (^s_) 

v r c,a ' 




in A 


(4l) 


RMS of Fit in A 


H/3 


]<3 - 


(1.666 - 0.324(Mg b)) 


X 


(0.838 + 0.076(Mg 6)) x 


0.351 


0.074 


H 7 A 


h = 


(1.666 - 0.324(Mg b)) 


X 


(0.847 + 0.149(Mg 6)) x 


0.165 


0.091 




Jj_ _ 

Fc,i 


(1.666 - 0.324(Mg b)) 


X 


(0.847 + 0.149(Mg b)) x 


0.165 


0.091 




Js_ — 

F c ,5 


(1.666 - 0.324(Mg b)) 


X 


(0.649 + 0.293(Mg 6)) x 


0.091 


0.139 


R5 F 


36 = 


(1.666 - 0.324(Mg b)) 


X 


(0.649 + 0.293(Mg b)) x 


0.091 


0.139 



Table 2: The range over which Mg b yields the most reliable results is 3.0 A < Mg b < 4.3 A . This 
Mg b range covers most elliptical galaxies. 



Figures 1 and 2 show the measurements of Virgo and SDSS galaxy averages along with the 
corrected SDSS data (light blue line). In all the graphs the Virgo data is in red while the SDSS 
data is in green. Also in Figures 1 and 2 model grids are plotted in blue and pink. The blue 
corresponds to models of various ages from 1.5 to 17 Gyrs and the pink corresponds to models of 
various metallicities from -2.00 to 0.50. For the H/3 correction the line fit is a little low for smaller 
galaxies with weaker Mg b, but not bad for the larger galaxies. See Figured! 
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Fig. 1. — Ha, H/5,and against Mg b indices for two data sets and models. The Virgo spectra 
(red symbols with error bars), the SDSS spectra (green symbols with error bars), single stellar 
population models (blue grid lines) from bottom to top of ages 17, 12, 8, 5, 3, 2, and 1.5 Gyrs and 
models from right to left of metallicites 0.5,0.25,0, —0.25, —0.5, —1, —1.5, and —2 (pink grid lines) 
are plotted. The fits to the index values for the SDSS galaxies after the correction for hydrogen 
emission are shown as a light blue lines in each panel. 

The fits for H7 and H<5 do not look as good when compared to the Virgo measurements (c.f. 
Figure 2). The reason is most likely the sensitivity of the H^y and H5 indices to changes in individual 
elemental abundances, with perhaps a non-negligible contribution from systematic errors between 
the two data sets. Elements N, C, and O have a profound and interdependent effect on this region 
of spectrum, and suggest themselves as candidates for further investigation. 
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Fig. 2. — H<5A,and H5p against Mg b indices for two data sets and models. The Virgo spectra 
(red symbols with error bars), the SDSS spectra (green symbols with error bars), single stellar 
population models (blue grid lines) from bottom to top of ages 17, 12, 8, 5, 3, 2, and 1.5 Gyrs and 
models from right to left of metallicites 0.5,0.25,0, —0.25, —0.5, —1, —1.5, and —2 (pink grid lines) 
are plotted. The fits to the index values for the SDSS galaxies after the correction for hydrogen 
emission are shown as a light blue lines in each panel. 



3. Discussion, Summary, and Conclusion 



Using the correction factor (j^-) as derived above and associatin g it with the linear offset 



in Ha this work has shown that for SDSS spectra cleaned like those in iGravesI (|2007l ). a emission 
correction factor on the order of 0.5 A in Ha and 0.2 A in H/3 needs to be applied to the quiescent 
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galaxies in order to better determine the mean age and metall icity of t hese g alactic averages. These 
correction factors are much larger than those estimated in I Graves! (|2007l ). where the estimated 
correction factors from measurements of Oil give Ha = 0.082 A and H/3 = 0.027 A. This work 
also shows that the higher-order Balmer lines H7 and H<5 have such modest emission corrections 
that other uncertainties dominate the error budget. At least for our comparison sample of Virgo 
cluster elliptical galaxies, H7 and H5 suffer from contamination from the varying of other elemental 
abundances, making age determination from these indices more complicated. This makes the need 
for correct Ha and H/3 measurements all the more important, since they are relatively insensitive 
to changes in metal abundances. 



These results are in qualitative agreement with the results from lEisenstein et al.l (|2003l ) who 
showed that in their SDSS spectra H/3 suffers from interstellar emission lines and speculated that 
non-solar abundance ratios were to blame for the differences in age determinations made from H/3, 
H7 and H<5. This is in agreement with these findings, where not only is there obvious hydrogen 
emission, but also it is clear that relative abundance ratios would certainly have to be taken into 
account in order to determine galactic ages from H7 and H5. 

For this work the Ha correction was chosen as the basis for the corrections due to the facts 
that Ha is 3 times more sensitive to hydrogen emission than H/? and the line fit for Ha as measured 
in the virgo galaxies had a much tighter fit RMS = 0.066 A compared with H/3 RMS = 0.15 A 
leading to a corrected SDSS line fit with much less uncertainty. Especially for galaxies inside the 
range 3.0 < Mg b < 4.3 A where the corrections for H/3 yield the most reliable results. Another 
reason for choosing Ha instead of directly using the correction that one could get from the H/3 plot 
is that we wanted to preserve H/3 for age determination. 

The plot for H/3 in figure 1 panel 2 still shows a discrepancy between the corrected SDSS data 
(light blue line) and the Virgo data (red line). This residual discrepancy could be due to a few 
variables such as varying abundance ratios within the virgo galaxies, slightly different decrements 
values than the ones used here, or there is the possibility that t he Ha emission corre ction may be 
in part due to an difference in the mean ages of the two samples. iThomas et al.l (|2005l ) showed that 
cluster galaxies tend to be around 2 Gyrs older than field galaxies. 

Another possibility is that Ha and H/? have different age sensitivities. In order to investigate 
the age sensitivities using the aforementioned model s the Z sensitivity parameter was calculated for 
both Ha and H/3 as it was in lWorthev et al.l (| 19941 ). The age sensitivity parameter is the ratio of 
the percentage change in age to the percentage change in Z of the index measured as shown below. 



_ [AI/(AZ/Z)] 

[AI/(Aage/age)} [ > 

Here AI /(AZ/Z) is the average of the changes in the measured index I as measured from the 
models where the age is 12 Gyrs and [Fe/H] changes form [Fe/H] = 0.25 to [Fe/H] =0.0 and [Fe/H] 
= 0.0 to [Fe/H] = -0.25 and Z= 0.0169 x 10^ H l Similarly, AI /(Aage/age) is the average of the 
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changes in the measured index I as measured from the models where the metallicity is solar and the 
age changes from 17 Gyrs to 12 Gyrs and 12 Gyrs to 8 Gyrs and age = 12 Gyrs. These sensitivities 
are shown in Table 3 along with the original IWorthev et al.l (| 19941 ) H/3 sensitivity. Note that the 
models indicate that both Ha and H/3 are age indicators of the same sensitivity. This means that 
the discrepancy between Ha and H/3 is more likely to be due to abundance ratios or the decrement 
values used than it is the age sensitivities of Ha and H/3. 



Table 3 



Index 


Z (sensitivity parameter) 


Ha 


0.8 


H/3 


0.8 


H/3 (Worthey et. al 94) 


0.6 



The applicability of this work for other grand SDSS averages may be limited due to the details 
of the sample selection. It is also unlikely to be applicable to individual red sequence galaxies due 
to the wide dispersion in index values and possible age e ffects. One possible application could be 



to scale those averages to match those of iGravesI (j2007l ). but it would be a better idea to echo 



the work shown here by comparing with the minimal-emission Virgo data set those averages and 
determining a new correction. 

It is conceptually possible to solve for both age and emission correction by considering both 
H/3 and Ha simultaneously, and increasing the emission correction until both indices give similar 
ages against a model grid. The obvious trouble with that is that the solution then becomes model 
dependent. There is also strong anticorrelation between derived age and emission correction. Fi- 
nally, there is the contribution of N emission near Ha, whose presence may cause spuriously large 
Ha emission measurements in individual galaxies. 

Speculation aside, it is safe to say that there is emission contamination in the SDSS spectra 
and that it is reasonably well accounted for by the linear fits presented in this paper. In data sets 
that include Ha, observed Oil and OIII do not need to be used as a proxy for Balmer emission. 
Emission contamination is much less of a problem for H7 and H<5, but interpretation of these indices 
is complicated by the probable effects of individual elemental abundances. 
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